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Abstract 

We consider the model of interaction between the immune system 
and tumor cells including a memory function that reflect the influence 
of the past states, to simulate the time needed by the latter to develop 
a chemical and cell mediated response to the presence of the tumor. 
The memory function is called delay kernel. The results are compared 
with those from other papers, concluding that the memory function 
introduces new instabilities in the system leading to an uncontrolable 
growth of the tumor. If the coefficient of the memory function is used 
as a bifurcation parameter, it is found that Hopf bifurcation occurs 
for kernel. The direction and stability of the bifurcating periodic so- 
lutions are determined. Some numerical simulations for justifying the 
theoretical analysis are also given. 
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1 Introduction 



As everyone knows, cancer is one of the most fearsome illness. It was 
declared the disease of 20th century. Many efforts were made to cure it, but 
to do this, first of all it is needed to understand its physiopathological me- 
chanisms. It was discovered that the human body is not completely helpless 
against this disease and it fights against cancer using its best and powerful 
weapon, namely immune system. 

In what follows, we will not make a general overview of the immune 
system, but we will mention briefly some of its components and aspects of 
the dynamics which appear in our model. The cell who performs directly 
the tumor elimination is T-lymphocyte, which is activated by b-lymphocytes 
through the cytokines [3,11]. It is not of less importance to mention the 
immuno depression, a phenomenon that appears in the tumor region, when 
the tumor increases its size, and leads to the deactivation of the lymphocytes. 

In the effort of modeling this process an important role plays the time 
delay. It is obvious for everyone that the biological process do not take 
place instantaneously and an amount of time is needed for, in our case, the 
interaction between immune system and the tumor [2,4]. During the years, 
some models, concerning tumor dynamics have been develop [6,7,12] and 
some of them includes time delay [1,10,13]. 

In what follows we propose a model of the interaction tumor-immune 
system using delay kernel. 

Let x(t) and y(t) denote respectively the number of malignant and lym- 
phocyte cells, for tel. The rate of malignant cells (x(t)) is given by [12]: 



We assume that the growth rate is proportional to x(t) and the decrease 
rate is proportional to the frequency of interaction with lymphocytes. The 
coefficients are a\ and a 2 , respectively, where a\ is tissue dependent. 
On the other hand, the growth rate of lymphocytes y(t) is described by [12]: 



It is proportional to the interaction with malignant cells and also to the flux 
per unit time of lymphocytes to the place of interaction. These effects are 
represented by the first and fourth terms in the right-hand side of equation 
(2). The mortality of the lymphocytes is proportional with y(t) (natural 



x(t) = a,ix(t) — a 2 x(t)y(t). 



(1) 



y(t) = hx(t)y(t) - b 2 x(t) - b 3 y(t) + b A . 



(2) 
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death) and also x(t), which express the immunodepression phenomenon. The 
term bix(t)y(t) is important for this study. It means the interaction between 
the two populations, x(t) and y(t) with a frequency b\ of recognition of 
malignant cells by the immune system. We consider the effect of influence 
of the past for this chemical signal mediated interaction which introduces 
the memory functions pi and p 2 , which are nonnegative bounded functions 
defined on [0, oo) and 



oo 

/ 



ki(s)ds = 1, / ski(s)ds < oo,i — 1,2. 



The evolution equations (1), (2) become now 
x(t) = aix(t) — a2x(t)y(t) 

oo oo 

y(t) = 6i( f k 1 (s)x(t - s)ds)( [ k 2 (s)y(t - s)ds) - b 2 x(t) - b 3 y(t) + 6 4 . (3) 



The memory functions are called delay kernels. The delay becomes a discrete 
one when the delay kernel is a delta function at a certain time. Usually, we 
employ the following form 

h(s) = igf + Ve-<^, <=1,2, 

for the memory function. When p = and p — 1, the memory functions are 

called "weak" and "strong" kernel, respectively. 

For ki(s) = 5(s — Ti),i = 1, 2, t\ > 0, r 2 > equation (4) is given by 

x(t) = aix(*) - a 2 x(t)y(t) 

y(t) = b x x(t - n)y(t - t 2 ) - b 2 x(t) - b 3 y(t) + b 4 . 

The model (4) with T\ = r 2 = r, is the model from [4] which has been studied 
using only numerical simulations. 

In this paper, we analyze the model (4) with the following initial values 
x 1 (9) = <p 1 (e), x 2 (9) = <p 2 (e), 9 e (-00,0} 
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and ipi , ipi as differentiable functions. 

The paper is organized as follows. In section 2, we discuss the local stabi- 
lity for the equilibrium states of system (4), for different forms of the delay 
kernels. We investigate the existence of the Hopf bifurcation with respect 
of the parameters of the delay kernels. In section 3, the direction of the 
Hopf bifurcation is analyzed by normal form theory and the center manifold 
theorem. Numerical simulations in order to justify the theoretical results are 
illustrated in section 4. Finally, some conclusions are made. 



2 Local stability and existence of the Hopf 
bifurcation 



We consider model (4) with parameters a±, a 2 , &i, b 2 , h, &4 assumed positives 

?>2 ^ b t ^ ai 
6l 63 (12 ' 

- (C\ 

f>3 



numbers and < || < ^. The equilibrium states of system (4) are the 
points L = (xq, yo) and L\ = (0, ¥), where 



b^ai — 6 4 a2 a\ 

aibi - a 2 b 2 a 2 

We analyzed the local stability in the equilibrium state L . We consider the 
following translation 

Xl (t) = x(t) - x Q , x 2 (t) = y(t) - y . (5) 

With respect to (5), the system (4) can be expressed as 

Xi(t) = -a 2 x x 2 (t) - a 2 xi(t)x 2 (t) 

00 00 
x 2 {t) = -b 2 xi{t) - b 3 x 2 (t) + bix J k 2 (s)x 2 (t - s)ds + biy Q J k 1 (s)x 1 (t - s)ds 



00 00 
+ 61 ( J h(s)x 1 (t - s)ds)( J k 2 (s)x 2 (t - s)ds). 


(6) 

The system (6) has 0=(0,0) as equilibrium state. 

To investigate the local stability of equilibrium state of the system (8), we 
linearize system (6). The linearized system of (6) is 

U(t) = AU(t) + B 1 U 1 (t) + B 2 U 2 (t), (7) 
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where 

*-UT>-(.I!).*-(S£) (8) 

with 

oo oo 

U(t) = ( Ul (t),u 2 (t)) T , Ui(t) = ( J k i (s)u 1 (t-s)ds, J k i (s)u 2 (t-s)ds) T , 1 = 1,2. 



The characteristic equation corresponding to system (7) is A(A)=0, where 

oo oo 

A(A) = det(XI ~ A -(J h(s)e- Xs ds)B 1 - ( J k 2 {s)e- Xs ds)B 2 ). (9) 



From (8) and (9), we have: 

oo oo 

A(A) = \ 2 + b 3 \-a 2 b 2 x + aib 1 x / k^e^^ds-Xb^o / k 2 (s)e~ Xs ds. (10) 



The equilibrium state Lq is locally asymptotically stable if and only if the 
eigenvalues of A (A) = have negative real parts. 

Because of the presence of two delay kernels k\ and k 2 in the equation A(A) = 
0, the analysis of the sign of real parts of eigenvalues is complicated and a 
direct approach cannot be considered. 

We analyze the eigenvalues for the equation A(A) = if the delay kernels k\ 
and k 2 are delta functions or k\ is delta function and k 2 is weak function. 
Using results from [3], we obtain: 

Proposition 2.1. // 

k^s) = 8{s - n), k 2 {s) = 5{s - r 2 ), n > 0, r 2 > (11) 

then 

(i) function (10) is given by 

A(A, n, t 2 ) = A 2 + b 3 X - a 2 b 2 x + a^x^-^ 1 - Xb lXo e~ XT2 ; (12) 
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(ii) if Ti = 0,r 2 = then the equilibrium state L of system (4) is locally 
asymptotic stable; 

(Hi) if 

0<t 1 + t 2 < , 

a\0\XQ 

then the equilibrium state L of the system (4) is asymptotically stable. 

Next, we study the existence of Hopf bifurcation of system (4) with k\ and 
k 2 given by (11), by choosing one of the delays as a bifurcation parameter, 
e.g. take T\ as the bifurcation parameter. First, we would like to know when 
A(A,Ti,t 2 ) = 0, where A(A, r 1 ,r 2 ) given by (12) has purely imaginary roots 
A = ±.iujQ{ujQ > 0) at Ti = tiq. Note that 

+ o, 2 b 2 XQ - a^xo cos(t 10 lu ) + biX L0 sin(r 2 cJo) = 
b 3 u - aiMo sin(ri u;o) - hx u cos(r 2 u ) = 0, 



(13) 



which implies that 

sin((no - r 2 )w ) = g(uo), 

where 

a/ - (b\xl -bl- 2a 2 b 2 x )uj 2 + [a\b\ - a\b\)xl 



2a 1 b\x\ 



to 



From (14), g'{uj) > 0. So g(u) is strictly monotonically increasing on [0, 00), 
with lim g(cu) = —00 and lim g(cu) = 00. Clearly, if r 10 > r 2 then g(u>) 

intersects sin((rio — t 2 )uj) only in a point. Hence A = iu>o is a simple root of 
equation A(A,Ti,t 2 ) = 0. Differentiating A(A,Ti,T 2 ) implicitly with respect 
to Ti, we obtain 

-^ e [^rJ A=jwo,n=no ~~ 

a{b\xl + aib\xlT 2 ujl cos((Tio-T 2 )cj )+ai&i^(u^ + a 2 6 2 £o) cos(ri cc;o) 

where 

k = 63 - aihxom cos(t 10 u ) - hx cos(t 2 uj ) + &i£ T2^o sin(r 2 cj ) 
l 2 = 2uj - ai6iX Ti sin(rioWo) + hx sin(r 2 cj ) + hx T 2 uj cos(r 2 a; ). 

From the above analysis and the standard Hopf bifurcation theory, we have 
the following result: 
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Proposition 2.2. Ifki(s) = S(s — Ti), k 2 (s) = S(s — r 2 ) and there is T\ = r 10 
for given r 2 > ; r 10 > r 2 so that equations (13) hold and 

. dX 

then a Hopf bifurcation occurs at L Q as T\ passes through r w . 
For given r 2 > 0, a solution for (13) is the pair (tio,u>io), where 

J\>7T 

rio = hr 2 , fc = 1,2, . . . 

and Wio is a positive root of the equation 

x 4 - (bjxl - &§ - 2a 2 & 2 x ):r 2 + (a 2 6 2 - a?6f )xq = 0. 
Proposition 2.3. If 

k^s) = 6(s-n), k 2 (s) = q 2 e~^ s ,r 1 > 0,g 2 > (15) 

then 

(i) function (10) is given by 

A(A,T!,g 2 ) = A 3 +p 2 A 2 + Pl A + p + (r 1 A + r )e- Ari , (16) 

u>aere 



p 2 = g 2 + & 3 , Pi = - a 2 b 2 x - bix q 2 
Po = -g2«2&2^o, n = a^ixo, r = aibix q 2 ] 



(17) 



(ii) if Ti =0 and 

4(ai&i - a 2 b 2 ) 2 < a 2 & 3 (f>i&4 - a 2 6 3 ) 

then for q 2 G (0, q , 2 i)U(a 22 , oo) the equilibrium state L is locally asymp- 
totic stable, where q 2 i,q 22 are the solutions of the equation 

(6 3 - b 1 x )x 2 + 63(63 - b l x )x + 6 3 (ai6i - a 2 b 2 )x = 0. 

Next, we study the existence of Hopf bifurcation for system (3) with p x 
and p 2 given by (15), by choosing the delay T\ as the bifurcation parameter. 
First, we would like to know when A(A, Ti, q 2 ) = 0, where A (A, Ti, q 2 ) is given 
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by (16), has purely imaginary roots A = ±zu; 01 (u;oi > 0) at T\ = Tu. Note 
that 



Po ~ P 2^01 + r o cos(a;oiTii) + r x sin(a; 01 rii) = 
- u\ x + pia; i + ri^oi cos^oiTn) - r sin(cj iTii) = 0, 



(18) 



which implies that 

< + (PS - 2pi)< + (pi - 2p p 2 + r?)< + p2 _ r 2 = o. (19) 

From (17), Pq < Tq and from (19) A = zc^oi is a simple root of the equation 
A(A,Ti,g 2 ) = 0. From (18) we obtain: 

1 ri^ i(P2^oi - Po) + ^(Pi^oi - o>oi) 
Tu = arctan — = — „— . 

uoi Pi^oi(^oi - Pi^oi) + r (Po - £>2^oi) 

Differentiating A(A, Ti,g 2 ) = implicitly with respect to r±, we obtain 

p ,_^A _ ^ i(riaWi - l 2 r ) 

ar\ m\ + m 2 

where 

™>i = (Pi ~ 3u; oi) cos(a;oirii) - 2p 2 u 01 sin(woirn) + r x - r T U 

m 2 = 2p 2 u 01 cos(a;oiTii) + (pi - 3^) sin(u;oiTii) - nm - rim^oi- 

Proposition 2.4. If k\(s) = 8(s — Ti), & 2 (s) = q 2 e~ q2S and T\ = Tu then 

dX 

Rc{~j^J \=iunn,ri=Tii T 1 

and a Hopf bifurcation occurs at L as T\ passes through T\\. 



3 Direction and stability of the Hopf bifurca- 
tion for ki(s) = S(s — ti), = £(s — T2) 

In what follows, we will study the direction and stability in two cases: in 
the first case the both kernels are delta function and in the second case the 
kernel k\ is delta function and the kernel k 2 is weak function. 

3.1. The case ki(s) = 8(s — T\), k 2 (s) = 8(s — r 2 ), n > 0, r 2 > 0. 
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In Proposition 12. II and 12.21 we obtained some conditions which guarantee 
that system (4) undergoes Hopf bifurcation at r = r 10 . In this section, 
we study the direction, the stability and the period of bifurcating periodic 
solutions. The method that we used is based on the normal form theory and 
the center manifold theorem introduced by [2]. 

From the previous section, we know that if r = r 10 , then all the roots of 
A(A, r 10 ,r 2 ) = 0, other than ±iu have negative real parts and any root 
of the form A(ti) = a{r{) + %uj{t\) satisfies a(r w ) = 0, uj(t w ) = u and 
da ^ ^ 7^ 0. For notational convenience let r% = t\q + fj,, /jL G R. Then \i = 
is the Hopf bifurcation value for (4). Without loss of generality, assume that 
t 10 > r 2 and define the space of C 1 functions as C 1 = C 1 ([— r 10 , 0], C 2 ). 
Suppose that for given a 1; a 2 , bi, b 2 , b 3 , 6 4 , r 2 , there is a 7*10 > at which (4) 
exhibits a Hopf bifurcation. In T\ = r 10 + fx, /x G M, we regard fi as the 
bifurcation parameter. For G C 1 , we define 




and 



r(h)4> = 



(o,of, ee[-r 10 ,o) 

(-02^(0)^(0), bi0i(-r lo )(/) 2 (-r 2 )f, = 



where 




A, # = 

b x 8(9 + t 2 ), ee[-r 2 ,o) 

B 2 S(9 + t w ), ee [-T W ,-T 2 ) 



and A, Bi, B 2 are given by (8). 

Then, we can rewrite (4) in the following vector form 



U t = A(fi)U t + RU t 



where 



U=( Ul ,u 2 ) T , U t = U(t + 6), 0G [-t 10) 0]. 
For tp G C 1 ([0, rio], C 2 ), the adjoint operator A* of A is defined as 
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dj>(s) 



, s e (0, no] 



For (/) e C([-no,0],C 2 ) and ip e C([0, r 10 ], C 2 ) we define the bilinear 
form 
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< >= V T (O)0(O) - J 1^(^-9)^(9)^)^, (20) 



-TIO ?=0 



where 77(0) = 77(0,0). 



Proposition 3.1. (%) The eigenvector of A(0) corresponding to eigenvalue 
Ai = 7cj groen fry 

M0) = K,f 2 ) T e Aie , 0e[-r lo ,O], 

7/j/iere 

a 2 b 2 - ai&ie A2T10 



ui = 1, v 2 



a 2 (h + Ai - 6iXoe A27 " 2 ) 
and A2 = Ai/ 

(m,) T/ie eigenvector of A* corresponding to eigenvalue \ 2 is 

h*(s) = ( Wl ,w 2 ) T e Xls , se[0,oo) 

where 



h 1 r 

Wi = — , W 2 = -, /1 = 



a 2 b 2 - ai6ie AlT1 ° 
77' ~' 77' a 2 Ai 
r/ = (fi + &i2/onoe AlT1() ) +^(1 + r 2 6 2 x e AlT2 ); 

(ra) H^i/i respect of (20), we have: 

<h*,h>=l, < h*,h >=< h*,h >= 0, <7f,7i>=l. 
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Using the approach of [2], [9] we next compute the coordinates of the 
center manifold Q at /j, = 0. Let X t = X(t + 9),9 e [— Tio, 0) be the solution 
of system (3) when /i — 0. 
Define 

z(t) =< h*,X t > w(t, 6) = X t - 2Re(z(t)h(6)). 
On the center manifold Q , we have 

2 2 

iw(t, 0) = w(z(t),z(t),9) = w 20 (6)j + Wll (9)zz + w 02 (9)j + ■■■ 

where z and z are the local coordinates of the center manifold Qq in the 
direction of h and h*, respectively. 
For the solution u t G VLq, we have: 



z(t) = \ 1 z(t)+g(z(t),z(t)), 

where 



z 2 z 2 z 2 z 



g(z, z) = g 20 — + guzz + g 2 i — + 92i^~ 

Proposition 3.2. For the system (4), the coefficients g 20 , gu, g Q2 , g 2 \ and the 
functions w 2 o(9), wu(9), wq 2 {9) are given by 



920 = W1/120 + ^2/220, 9u = + w 2 f 211 , 

902 = W1/102 + W2/202, #21 = W1/121 + w 2 f 221 . 



(21) 



where 



120' 



220 



/120 = -2a 2 viw 2 , /in = -2a 2 Re(viv 2 ) , f W2 = f 

/220 = 26 lVl t;2e A2{Tl0+T2) , /211 = 26 1 i?e(tw3 Air2+A2 ^) ) /202 = 7 
/121 = -02(2^1^211(0) + W[w 220 (0) + 2v 2 w ul (0) +u[w 120 (0)) 
/221 = &i(2t;ie A2rio w; 211 (-r2) +^e AlT10 «;22o(-r2) +2 V2 e A2r2 wiii(-r 10 ) + 
+ ^e Air2 «; 1 2o(-Tio)) 
and 

W2o(0) = (w 120 (9),w 220 (9)) T , w n (9) = (w U i(9) , w 211 (9)) T , 
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and 
where 



E 
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Ai 3Ai 

Wn (6) = ^h(0)e^ 9 - 9 -^Wy X2d + E 2 , 
Ai Ai 



Ei = (E n , E 12 ) T , E 2 = (E 2 i, E 2 2) T , 

(2Ai + 6 3 - bix e AlT2 )/i2o - 02^0/220 

2Ai(-2Ai - b 3 + hxoe 2 ^) + a 2 x (-b 2 + hy e 2X ^) 



_ 2XiE n + /120 rp _ (bix — b 3 )E 22 + /211 /in 

t^\2 — j t^2\ — ; '-i j -^22 — ■ 

a 2 x ~ b 1 y -b 2 " a 2 x 

Based on the above analysis and calculation, we can see that each g^ in 
(21) is determined by the parameters and delays in system (3). Thus, we can 
explicitly compute the following quantities: 

. . % . . .0 1 . .9 X i?21 

Wo(0) = 2^(#20#n ~ 2 l^iil ~ 3I502I ) + -y 
^eCio(O) 

^ 20 " iteA'(O) (22) 
/mC 10 (0) + /i 20 ImA'(0) 



T20 — 

^0 

/3 20 = 2Re (Cio(O)). 

In summary, this leads to the following result: 

Theorem 3.1. In formulas (22), fi 2 o determines the directions of the Hopf 
bifurcations: if ' fi 2 o > 0(< 0) the Hopf bifurcation is supercritical (subcritical) 
and the bifurcating periodic solutions exist for n > Tio(< no)/ ^20 determines 
the stability of the bifurcation periodic solutions: the solutions are orbitally 
stable (unstable) if /3 2 o < 0(> 0) and T 20 determines the periodic solutions: 
the period increases (decreases) ifT 2 o > 0(< 0). 
In (22) Re(X'(0)) and Im(X'(0)) are given by 

Re(X'(0)) = R e (-^-) x=i ^ T=T10 
dX 

Im(X'(0)) = Im(—) x=iui0jT=T10 
12 



where 

d\ ai6iXoAe~ An 

dn b 3 + 2\- ai&iX Tie- Ar i - hx (l - Ar 2 )e- Ar2 ' 

3.2. The case fci(s) = 8(s - n), fc 2 (s) = g 2 e~ 92S , ri > 0, g 2 > 0. 

For fci(s) = S(s — Ti), k 2 {s) = q 2 e~ q2S , T\ > 0, q 2 > 0, system (6) is 
given by: 



xi(t) = -a 2 x x 2 (t) - a 2 x 1 (t)x 2 (t), 

x 2 {t) = -b2Xi(t)-b 3 X2(t)+bix Q x s (t)+biyoX2(t-T 1 )+bix s (t)x 2 (t-T 1 ), 
X3(t) = q2{x 2 (t) - x 3 (t)). 



We linearize system (23) and obtain: 

V(t) = A 1 V(t) + dV(t - n) 



(23) 



where 





(° 




-a 2 x 







-b 


2 


-h 


bix 




\o 






-Q2 






(0 







Ci 







hyo 









1° 





V 



with V(t) = ( Ul (t),u 2 (t),u 3 (t)) T . 

The characteristic equation of system (23) is given by A(\,ri, q 2 ) = 0, 
where A(X,Ti,q 2 ) is function (16). We consider T\ = T\\ + /i, ji G IR and 
C 1 = C 1 ([— Tn, 0], C 2 ). We regard /x as the bifurcation parameter. Then, for 
4> E C 1 , we define 



Ai(fj,)<f> = 



d<t>(e) 
de 



, fe[-ru,o) 



/ drj(t,^)4>(t), 6 = 



-Til 



and 



13 



R(llU -S (0,0, of, ee[-m,0) 

lW "1 (-a 2 0i(O)0 2 (O),6 1 03(O)02(-r 11 ),Of, # = 

where 

J A, 6» = 

\ d^ + m), ^[-t U) 0). 

Then, we can rewrite (23) in the following vector form 

U t =A 1 ( fJl )U t + R 1 U t , 

where 

U t = U(t + 8), 0e[-rn,O]. 
For -0 G C 1 ([0, Tn], C 2 ), the adjoint operator .A^ of .Ai is defined as 



dip(s) 
ds 





, s e (o, m] 



/ drj T (t,0)i(;(-t), s = 0. 

-Til 

For (p G C([-Tn,0],C 2 ) and ^ e C([0,r u ],C 2 ) we define the bilinear 
form 

o e 

<0,^>=^(O)0(O)-| J^ T (^-9)dr ] (e)<j ) (0dt (24) 



-rii £=0 



where rj(9) = 7/(0, 0). 



Proposition 3.3. (i) The eigenvector of Ai(0) corresponding to eigen- 
value Ai = iujQi is given by 

h{9) = (v 1 ,v 2 ,v 3 ) T e Md , 9 G [—Tii, 0], 

where 

vi = (Ai + g 2 )(Ai + 6 3 -^i2/oe A2T11 )-g2&i^o, v 2 = -b 2 (X 1 + q 2 ),v 3 = -b 2 q 2 
and X 2 = Ai; 
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(ii) The eigenvector of A\ corresponding to eigenvalue X 2 is 
h*(s) = {w 1 ,w 2 ,w 3 ) T e XlS , se[0,oo) 

where 

fi 1 h t b<2 t blX ° 

Wi = — ,W 2 = ~,W 3 = —,f 1 = - — ,/ 3 



V V V A 2 A 2 + q 2 

V = fivi + Ml ~ - e Aim - A 2 rn6 2 e Airi1 )) + f 3 W 3] 

X 2 

(Hi) With respect to (24), we have: 

<h*,h>=l, < h*,h >=< 7T, h >= 0, <h*,h>=l. 

Using the approach of [2], [9] we next compute the coordinates of the 
center manifold f2 at /i = 0. Let X t = X(t + 9),9 G [— Tu, 0) be the solution 
of system (3) when \x = 0. 
Define 

z(t) =< h*,X t > w(t, 6) = X t - 2Re(z(t)h(9)). 
On the center manifold flo, we have 

2 2 

w(t, 0) = iu(z(*), S(t), 0) = w 20 (^)y + + w m {6) Z - + ■■■ 

where z and z are the local coordinates of the center manifold Qq in the 
direction of h and h*, respectively. 
For the solution X t £ Q , we have: 



z(t) = X 1 z(t)+g(z(t),z(t)), 

where 

-j^ z^~z 
g(z,z) = 920— + Qiizz + g 2 i— + g 2 i — 

Proposition 3.4. For the system (23), the coefficients g 2 Q, gu, go 2 , g 2 \ and 
the functions w 20 (9) , Wu(9) , w 02 (9) are given by 

920 = wi/120 + ^2/220 + w 3 f 320 , gu = wi/in + W2/211 + W3/311, ^ 
902 = wi/102 + ^2/202 + w 3 f 302 , 921 = wi/121 + W2/221 + w 3 f 321 . 



15 



where 

/120 = -2a 2 viv 2 , /m = -2a 2 Re(v 1 vi), f W2 = f 120 , 

/220 = 2b lV2 v 3 e x ^\ f 2U = 2b 1 Re(v 2 v 3 e x ^),f 202 = J 220 

/320 = /311 — /302 = 



/121 = -02(2^1^211(0) + ^1^220(0) + 2v 2 w ul (0) + v 2 w 120 (0)) 

/221 = &i(2u 3 «;2ii(-Tn) +v^w 220 (-Tn) + 2v 2 e X2Tll w 3n (0) + v^e XlTll w 320 (0)) 

/321 = 0; 

and 

w 2Q (9) = (w 120 (6),w 220 (6),w 320 (6)) T , w u (0) = (w lu (6), w 211 (6), w 3U (8)) T , 



and 

where 
En 

Ei 3 



W20 (9) = -^h(0)e x ^ - I^Wy* 9 + Eie^° 

Al OAl 

Wll (9) = ^h(0)e Me - ^-h{0)e X2e + E 2 , 
Ai Ai 

Ei = (En, Ei 2 , Ei 3 ) T ', E 2 = (E 2 i, E 22 , E 23 ) T , 

a 2 x§Ei 2 + fi 20 _ 2A1/220 — b 2 fi 2 Q 

2A^ ' 12 ~ 2Ai(2Ai + 63 - 6n/oe 2A ^") - a 2 W 
Q2E11 (63 - biy )E 22 - /211 /m 

7TT ; ; ^21 — } 1 -&22 — , -&23 — ^21- 

2Ai + q 2 b 2 a 2 x 

We can explicitly compute the following quantities Cn(0), /i2i, T 2 i, /32i : 

Cn(0) = ^^(^20^11 - 2|^n| 2 - ^ I ^02 1 2 ) + ^ 
ReCn(0) 

^ ~ ~ ReX'(0) (26) 

/mCii(0)+^2i/mA'(0) 
1 21 — 

/3 21 = 2Re(C u (0)). 
In summary, this leads to the following result: 
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Theorem 3.2. In formulas (26), fi 2 i determines the directions of the Hopf 
bifurcations: if > 0(< 0) the Hopf bifurcation is supercritical (subcritical) 
and the bifurcating periodic solutions exist for n > ni(< m)/ (3 2 i determines 
the stability of the bifurcation periodic solutions: the solutions are orbitally 
stable (unstable) if fi 2 \ < 0(> 0) and T 2i determines the periodic solutions: 
the period increases (decreases) ifT 21 > 0(< 0). 
In (26) Re(X'(0)) and Im(X'(0)) are given by 

dX 

Re(X'(0)) = Re(—) x=iU0UT=ni 
Im(X'(0)) = Im(^-) x=kJ01:T=T11 

where 

dX (nA 2 + r X - ri)e~ Ari 



dri 3A 2 + 2p 2 X + p 1 — (r x A + r )ri 



4. Numerical simulations. 

For the numerical simulations we use Maple 9.5. In this section, we 
consider system (6) with a\ = 2.5, a 2 — 1, b\ — 1, b 2 = 0.4, 63 = 0.95, 64 = 2. 
We obtain: x = 0.1524390244, y = 2.5. 

In the first case, ki(s) = 5(s — Ti), k 2 (s) = S(s — t 2 ), for r 2 = 0.01, 
we have: u = 0.6124295863, fi 2 = 630.5712553, (3 2 = 125.5070607, T 2 = 
10.25944116, no = 9.541873607. Then the Hopf bifurcation is supercritical 
and the bifurcating periodic solutions exist for r > r 10 ; the solutions are 
orbitally unstable and the period of the solution increases. The waveforms 
are displayed in Figl and Fig2 and the phase plane diagram of the state 
variables x(t), y(t) is displayed in Fig3: 



17 



Fig.1.(t,x(t» 



Fig.2.(t,y(t» 




In the second case, ki(s) = 5(s — Ti), /^(s) = qie Q2S for q 2 = 0.1, we 
have: u; i = 0.2235621332, = 7.926079992, /3 21 = 0.04097046568, T 21 = 
0.3275619874, = 10.38589492. Then the Hopf bifurcation is supercritical 
and the bifurcating periodic solutions exist for T\ > r* : ; the solutions are 
orbitally unstable and the period of the solution increases. The waveforms 
are displayed in Fig5 and Fig6 and the phase plane diagram of the state 
variables x(t), y(t) is displayed in Fig7: 
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Fig.4.(t,x(t)) 



Fig.5.(t,y(t)) 




0.1515 0.152 0.1525 0.153 0.1535 



For q 2 = 0.1, we have: uj 01 = 0.9506753825, fi 21 = -0.6058263333, foi = 
-0.001118156944, T 21 = -0.07864963978, r n = 23.03933807. Then the Hopf 
bifurcation is subcritical and the bifurcating periodic solutions exist for T\ > 
Tu] the solutions are orbitally stable and the period of the solution decreases. 
The waveforms are displayed in Fig5 and Fig6 and the phase plane diagram 
of the state variables x(t), y(t) is displayed in Fig7 and Fig8: 
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Fig.7.(t,x(t)) 



Fig.8.(t,y(t)) 




0.1515 0.152 0.1525 0.153 0.1535 



5. Conclusions. 

This paper was focused on mathematical analysis of a model which describes 
the interaction between immune system and the tumor cells. The model is 
an improved one by using the delay kernel. Taking the average time delay as 
a parameter, it has been proved that the Hopf bifurcation occurs when this 
parameter passes through a critical value. In a future work it will be studied 
the mathematical aspects of the effect of immunotherapy on the development 
of the cancer. 



20 



References 



[1] N. Buric, D. Todorovic, Dynamics of delay- differential equations model- 
ing immunology of tumor growth, Chaos, Solitons and Fractals 13(2002), 
645-655. 

[2] H.M. Byrne, The effect of time delays on the dynamic of avascular tumor 
growth, Math. Biosciences 144(1997), 83-117. 

[3] S. Dezfouli, I. Hatzinisiriou, S.D.J. Ralph, Use of cytokines in cancer 
vaccines /immunotherapy: recent developments improve survival rates 
for patients with metstatic malignancy, Curr. Pharm. Dcs., 11(2005), 
3511-3530. 

[4] M. Galach Dynamics of the tumor-immune system competition- the ef- 
fect of time delay, Int. J. Appl., Math. Comput. ScL, 13(2003), 395-406. 

[5] B.D. Hassard, N.D. Kazarinoff, H.Y. Wan, Theory and Applications of 
Hopf Bifurcation, Cambrige University Press, Cambrige (1981). 

[6] D. Kirschner, J.C. Panetta, Modelling imunotherapy of the tumor- 
immune system interaction, J. Math. Biol., 38(1998), 235-252. 

[7] V.A. Kuznetsov, I. Makalkin, M.A. Taylor, A.S. Perelson, Nonlinear 
dynamics of immunogenic tumors: parameter estimation and global bi- 
furcation analysis, Bull. Math. Biology, 56(1994), 295-321. 

[8] X. Liao, G. Chen, Hopf bifurcation and chaos analysis of Chen's system 
with distributed delays, Chaos, Solitons and Fractals, 25(2005), 197-220. 

[9] G. Mircea, M. Neam^u, D. Opri§, Dynamical systems from economy, 
mechanic and biology described by differential equations with time delay, 
Editura Mirton, 2003 (in Romanian). 

[10] D. Rodriguez-Perez, O. Sotolongo-Grau, R. Espinosa Riquelme, O. 
Sotolongo-Costa, J. A. Santos Miranda, J.C.Antoranz, Tumors un- 
der periodic theraphy- Role of the immune response time delay, 
|arXiv:q-bio.fO /0604008 vl 7 Apr 2006; 

[11] F. Saleh, W. Renno, I. Klepacek, G. Asfar, H. Dashti, P. Romero, 
A. Dashti Melanoma Immunotherapy :past, present and future, Curr. 
Pharm. Des., 11(2005), 3461-3473. 



21 



[12] O. Sotolongo-Costa, L. Morales Molina, D. Rodriguez Perez, J. C. An- 
toranz, M. Chacon Reyes, Behavior of tumors under nonstationary the- 
raphy, |a7Xw:math-ph/0203057| v5 26 Apr 2002. 

[13] M. Villasana, A. Radunskaya, A delay differential equation model for 
tumor growth, J. Math. Biol. 47 (2003), 270-294; 



22 



